- 1—1 

X 



|hep-ph/ 0006322 

YUMS 00-05 



Hadronic Invariant Mass Spectrum 
in B — > X u lv Decay with Lepton Energy Cut 

O ' 
O 
o , 

CN , K. K. Jeong a , C. S. Kim a ' b * and Yeong Gyun Kim 

G . 

a : Department of Physics, Yonsei University and IPAP, Seoul 120-749, Korea 

l—i '. 

00 ' b : Department of Physics, University of Wisconsin, Madison, WI 53706, USA 

(N ' 

c : YITP, Kyoto University, Kyoto 606-8502, Japan 
> ! (February 1, 2008) 

(N 

m 
\o 
o 
o 
o 

a: 



Abstract 



We discuss the implications of charged lepton energy cut to the hadronic in- 
variant mass spectrum in charmless semileptonic B decays. Charged-lepton 
energy cut is inevitable in order to remove secondary leptonic events such as 



5— i 

■ b — > c, r — > I, and to identify the chaged leptons at detectors experimentally. We 

consider three possible lepton energy cuts, Ef nts = 0.6,1.5,2.3 GeV, and found 
that with the most probable cuts Ef nt = 1.5 GeV and M^ ax = 1.5 (1.86) GeV, 
45 ~ 60% (58 ~ 67%) of decay events survive. Therefore, B — ► X u lu decay 
events can be efficiently distinguished from B — > X c lv decay events. We also 
discuss the possible model dependence on the results. 
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I. INTRODUCTION 



The determination of the CKM parameter \V u b\ is important for constructing the so-called 
unitary triangle. It is hard to determine \V u b\ from semileptonic B-meson decays because 
Cabbibo dominated decay mode B —>■ XJu obscures the B —>■ X u lv mode. The traditional 
method for extracting |V^&| from experimental data involves a study of the charged lepton 
energy spectrum in inclusive semileptonic B decays, B — ► X u c lv Q. The b —>■ u events are 
selected above charm threshold, i.e. for lepton energies E\ above (Mj, — M D ) / (2Mb) ~ 2.3 
GeV. However, this cut on Ei is not very efficient (only less than 10% of b — > u events survive) 
and also the dependence of the lepton energy spectrum on perturbative and non-perturbative 
QCD corrections is the strongest in this end-point region 

As an alternative, the determination of \V u b\ may come from the measurement of the 
hadronic invariant mass spectrum [0] in the region Mx < Md. For B — ► X c lv decays, one 
necessarily has Mx > Md = 1.86 GeV. Therefore, if we impose a condition Mx < Md, the re- 
sulting events come only from B — > X u lv decay and most of the B — > XJlv decays are expected 
to lie in this region. There is an experimental problem to be expected, though - the charm- 
leaking of misidentified charmed particles below the kinematic b — *> c threshold. To avoid this 
leakage, one may concentrate on hadronic invariant mass below a certain value Mj£ ax (< Md), 
say Mx ax = 1.5 GeV. The detailed studies of this method were performed in Refs. [[SHlOU- The 
integrated fraction of events was introduced, 



and studied in its sensitivity to the three basic parameters, /i^, nib and a s . 

However, we cannot apply the above results to the experimental data directly. From the 
practical point of view, the leptons with low energy, i.e. less than 0.6 GeV, cannot be experi- 
mentally identified within the dectectors. And also a larger lepton energy cut might be needed 
to select 'prompt leptonic events' (6 — > I) from 'secondary leptonic events' (b — >• c — > I, r — > I). 
An experimental method such as the technique of neutrino reconstruction |11[| , which would be 
used to measure hadronic invariant mass directly and inclusively, may require a lower cut on 
the charged lepton energy. Therefore, the hadronic invariant mass spectrum would be affected 
by the various lepton energy cuts. In this Letter, we study the effects of lepton energy cut on 




(1) 



the hadronic invariant mass spectrum in inclusive charmless semileptonic B decays and discuss 
their implications. 



II. DIFFERENTIAL DECAY RATE 

At the parton level, the most general hadronic tensor for B — > X u lv is given by 

Wn„(p, v) =Wi(v- p, p 2 ) {p^Vy + p u v^ - g^v ■ p - it^ a pp a v p ) 
- W 2 (v -p,p 2 )g^ + W 3 (v ■ p,p 2 )Vfj,v u 

+ W 4 (v -p,p 2 )(p^ +p u v ll ) + W 5 (v -p^PyPv . (2) 

Here v is the 6-quark velocity and p is the total parton momentum. The total momentum 
carried by the leptons is q = m^v — p. At the tree level, W\ = 25(p 2 ) and all other = 0. 
We follow the 0(a s ) corrections to the hadronic tensor from the paper of De Fazio and Neubert 



1C]. Introducing the scaling variables 



2#i ~ 2 P 2 A 2v-p 

x = , p = — x and z = , Id) 

rrib mf nib 

where E\ is the charged-lepton energy which is defined in the S-meson rest frame, the triple 
differential decay rate is given by 

, = 12r { (l +x -z)(z-x- f^W, + (l-z + f)^W 2 

dxdzap z I 2 2 

+ (x(z -x)- P^iWs + 2m b W i + mlW 5 )}, (4) 

where x = 1 — x, and 

G 2 F \V uh \ 2 ml 

10 = 192*3 • (5) 
In terms of the parton variables, the total invariant mass of the hadronic final state is given 

by 

s H = p 2 + 2Av ■ p + A 2 , (6) 
where A = Mb — m^. Using relation (6), and the scaling variables, 

s H = s H /ml and e = A/m^ ; 



we find the following double differential decay rate 



d 2 r r , d 3 r 

dz- 



p = s H — ez — e 



dx ds H J dx dz dp 2 
where the phase space for the relevant variables is given by 



(7) 



e + x ~~ ~~ e 
efl + e)-s 



- < x < 1, (8) 



for e 2 < s H < e(l 



■ 



< z < 1 - e + 



e + x 1 + e 

0< x< (1 + £)2 -^ , (9) 
for e(l + e) < s H < (1 + ef . 

If we integrate the above double differential decay rate over the variable x, we get the single 
differential decay rate for s H 

dr ^^L, do) 



ds H J dxds H 



which reproduces numerically the result of Ref. [Kj . 

In order to obtain the physical decay distributions, we should also consider the non- 
perturbative corrections. The physical decay distributions are obtained from convolution of 
the above parton level spectra with non-perturbative shape function F(k + ), which governs the 
light-cone momentum distribution of the heavy quark inside the i?-meson ||||. The convolu- 
tion of parton spectra with this function is such that for the decay distributions the 6-quark 
mass m;, is replaced by the momentum dependent mass, m& + k + , and similarly the parameter 
A = Mb — rrib is replaced by A — k + . Here k + can take values between -mj and A, with a 
distribution centered around k + = and with a characteristic width of 0(A). Then, the scaling 
variables x, s H and e are replaced by the new variables 

where q + = A — k + . The physical double differential decay rate for the charged- lepton energy 
and the hadronic invariant mass is given by 



d 2 r KT, F(A-g + ) dr A 

" 2 /_ d ^+— — z^rK'^ 6 ?)- ( 12 ) 



d^ds H jo (M B — q + ) 3 dxds H 

where g™ ax = min(Ms — 2E[, y/s^). Finally, the physical distribution for the hadronic invariant 
mass can be obtained in two different ways; 



dr /• d 2 r 



d.s„ J dEi ds H ^ 

, F(A-q+) dr , A . 
o (M B - q + y ds H 

which should give the same results. We note that, after the implementation of such Fermi 
motion, now the kinematic variables take values in the entire phase space determined by hadron 
kinematics, i.e. < E x < M B /2 and < s H (= M|) < M 2 B . Eqs. (12-14) are our starting 
point for the numerical calculations. 



III. NUMERICAL ANALYSES AND CONCLUSIONS 



To perform the numerical calculation we should choose a specific form for the shape function 
F(k+ = A — q + ). It is subject to the constraints on the moments of the function 



A n = (k») = J dk + klF(k A 



which are given by the expectation values of local heavy quark operators. In practice we know 
only the size of the first few moments; one finds 

1 o 



A = 1, Ai = Q and A 2 = -/x; 



(15) 



where is the average momentum squared of the b quark inside the 73-meson. The parameters 
A and were obtained by the HQET and QCD sum rules in Refs. |12|,l3|; 



A = 0.4 ~ 0.6, 
/i 2 =0.6±0.1 [12] or /i 2 = 0.10 ±0.05 M. 



(16) 



One then chooses some reasonable ansatz for F(k + ); its parameters are adjusted so as to 
reproduce its known moments. Several functional forms for this function have been suggested 
in the literature |6l,|T4T-[T6|| • We adopt the simple form of Ref . |l6j , 



F(k+) = N(l - x) a e {l+a)x with x = 5t < 1, (17) 

which is such that A% = by construction (neglecting exponentially small terms in m^/A), 
whereas the condition for A Q fixes the normalization N. The parameter a can be related to 
the second moment, yielding A2 = A 2 / (1 + a). A typical choice of values is rrib = 4.8 GeV and 
a = 1.29, corresponding to A « 0.48 GeV and /x 2 ~ 0.3 GeV 2 . Below, we keep the parameter 
'a' fixed, and consider the three choices m b = 4.65, 4.8 and 4.95 GeV. These choices correspond 
to the following sets, 

(A(GeV),/i 2 (GeV 2 )) = (0.63,0.52), (0.48,0.30) and (0.33,0.14), (18) 

respectively. We fix the QCD coupling a s = 0.22. 

Now the calculation of the hadronic invariant mass spectrum is straightforward. The 3- 
dimensional plot of double differential decay rate, 

^ 2M X -^- (19) 



dM x dEi ds H dE l 

is shown in Fig. 1, with mb = 4.8 GeV. The lines in the bottom plane show the contour lines. 
The three straight lines in contour figure, which are the projections of three curved lines in 
the 3-dimensional figure, indicate the charged-lepton energy cuts, i.e. Ef nt = 0.6, 1.5, 2.3 GeV, 
respectively. From the figure, one can find that a significant fraction (~ 90%) of the decay 
events lie below the charmed states, i.e. have Mx < 1-86 GeV. This is in sharp contrast with 
the usual cut on charged-lepton energy E\. Due to measurement errors, there will be a tail 
from b —>■ c decays below Mp. To avoid this leakage, one can concentrate on hadronic invariant 
masses below a certain value M^ ax < M D . With a lower cutoff on M x (say Mj£ ax = 1.5 GeV 
or 1.6 GeV), the majority of decays still appears below the M^ ax region. 

But in reality we cannot apply the above results to experimental data directly. In B — > Xlv 
semileptonic decay experiment, the secondary electrons such as b — > c — > I cascade decay 
contaminates the signal of the pure B — > XJu decay. The secondary electrons have typically 
lower energy than the primary electrons because the mass of D-meson is much less than the 
mass of .B-meson. Therefore, only the electrons with energy greater than 1.5 GeV can be 
considered to be purely from B — > X cu lv. Using the double lepton tagging method, one can 
reconstruct the primary leptons with energy less than 1.5 GeV |fL7| . But in this method, both 



5-mesons from T(45 l ) should decay semileptonically. So it gives much less data samples. Also 
the experimental method such as neutrino reconstruction technique, which would be used to 
measure the hadronic invariant mass inclusively, requires a lower cut on the charged lepton 
energy. Therefore, we apply the charged-lepton energy cuts Ef ut = 0.6,1.5,2.3 GeV on the 
hadronic invariant mass spectrum. The cut Ef nt = 0.6 GeV corresponds to the double lepton 
tagging data, and Ef ut = 1.5 GeV corresponds to the normal single lepton tagging data. We 
also apply Ef ut = 2.3 GeV to compare the data with the previously measured charged lepton 
energy data as the kinematic boundary to separate the signal from b — > civ decay. 

Fig. 2 shows the hadronic invariant mass distribution with varying nib mass as 4.65, 4.8, 
4.95 GeV. It is normalized in units of (|T4&| 2 x GeV). Note that in the b — > u semileptonic decay, 
the leptons have quite large energies so that M x distribution with Ef nt = 0.6 GeV is not much 
different from the distribution without any lepton energy cut. In the case of Ef ut = 2.3 GeV, 
we can see that there is no decay event above M x > 1.86 GeV. This is consistent with the 
kinematics that leptons with energy greater than 2.3 GeV are purely from b — > ulv decay. And 
also we can see that this lepton energy cut is very inefficient because only a small part of the 
decays survive this cut, as mentioned earlier. 

Fig. 3 shows the integrated decay rate up to Mjp", 



The numerical values are summarized in Tables I, II and III. Here we can see that if we 
consider total decay rate with Ef nt = 2.3 GeV, only 9 ~ 14% of decay events remain. When 
we use M™ ax = Md cut, 82 ~ 94% of decay events remain without any lepton energy cut, and 
79 ~ 92% (58 ~ 67%) of decay events survive for Ef nt = 0.6 (1.5) GeV. In this case, there might 
be a difficulty to draw the reliable conclusion because of the possible charm-leaking. However, 
even in this case, the hard lepton energy cut may be helpful to reduce the contamination from 
the charm-leaking because the lepton energy of B — > XJu decays is softer compared to that 
of B — > X u lv decays. If we further reduce the maximum value of Mx, i.e. Mjp 1 * = 1.5 GeV, 
60 ~ 81% of the events remain without lepton energy cut, and 59 ~ 80% (45% ~ 60%) 
of the events survive for Ef ut = 0.6 (1.5) GeV. As we mentioned previously, in the case of 
Ef ut = 0.6 GeV, both _B-mesons should decay semileptonically to be identified as prompt 
leptonic events. Since Ef ut = 2.3 GeV gives only ~ 10% of the total decay rate, the cut 




(20) 



Ef nt = 1.5 GeV would be our best option. 

Now we discuss the dependence of the results on the various input parameters and the 
choice of the universal distribution function. The dependence on /x 2 is found not so significant, 
compared to the A (or equivalently m&) dependence. The decay rates with the parameters 
m& = 4.8 GeV and /i^ = 0.6 GeV 2 are shown in Table IV. By comparing Table IV with Tables 
I, II and III, we can see that the main uncertainty in decay rates comes from the uncertainty 
in A (or m^). The dependence of the result on the a s variation could be quite large since the 
perturbative correction to the total decay rate is linear with the parameter a s , and the size of 
a s correction is almost 20% of the leading approximation. From Tables I, II, III and V, we 
can see that the effects of the variation of a s from 0.22 to 0.35, with fixed m& = 4.8 GeV, are 
similar to the case of the variation of rrib from 4.8 to 4.65 GeV, with fixed a s = 0.22. 

Next, in order to estimate the possible dependence of the results on the choice of universal 



distribution funtion, we adopt ACCMM [18| model- induced distribution function [|6|,|T5l 

1 1 



ACCMM 
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_ Q± 






v F q+ 


Pf. 





(21) 



As seen, this shape function is dependent on the two parameters, p F and m sp . We choose the 
three sets, 

(p F ,m sp ) = (0.504,0.251), (0.383,0.193) and (0.262,0.136) (in units of GeV), 

which respectively correspond to the above three choices of mj = 4.65, 4.8 and 4.95 GeV with 
fixed 'a' of Eq. (17). The integrated decay rates T(M™ ax ) are summarized in Tables VI and 
VII. The resulting decay rates are almost the same as the corresponding results of Tables I and 
III. The differences in the decay rate for two universal distribution functions are within 2% 
only. 

Finally we note that the results would be affected by some resonance effects around M™ ax . 
The real result would be a sum of all the exclusive decays in which a few resonances dominate 
at some specific Mx's, so the actual M x distribution will be with humps and bumps, while our 
results are smoothed inclusive results using the duality. However, our integrated results would 
be quite correct, once there is no significant resonance around the region of Mx cut. If M]f ax is 
around 1 GeV, then there are a few important resonances, but if Mx is large enough, e.g. 1.5 



GeV, then there is no significant resonance from b — > u, but rather there will be decays with 
many pions. 

In summary, we investigated the effects of charged lepton energy cut on the hadronic in- 
variant mass spectrum for B — > X u lu decays and their implications. As is well known, the 
charged-lepton energy cut is experimentally inevitable in order to remove secondary leptonic 
events such as b — > c — > I, r — > I, and to identify the chaged leptons at detectors. We found 
that with Ef ut = 1.5 GeV and M™ ax = 1.5 (1.86) GeV, 45 ~ 60% (58 ~ 67%) of decay events 
survive, and with Ef nt = 0.6 GeV and Mj? ax = 1.5 (1.86) GeV, 59 ~ 80% (79 ~ 92%) survive. 
Finally without Ef ut , 60 ~ 81% (82 ~ 94%) of decay events survive with Mf ax = 1.5 (1.86) 
GeV. Therefore, B — > X u lv decay events can be efficiently distinguished from B — > X c lu decay 
events by using the hadronic recoil mass and the charged lepton energy together. 
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TABLES 

TABLE I. 

Ltotai with charged-lepton energy cuts with fixed a s — 0.22 and a — 1.29. The values 
of raj, = 4.65, 4.8, 4.95 GeV correspond to the set of values of Eq. (18) for A, 

mft(GeV) Ef^iGeV) L total (|\/J 2 GeV) percentage 



4.8 


o 


0.503 


100% 




0.6 


0.489 


97% 




1.5 


0.336 


67% 




2.3 


0.054 


11% 


4.65 





0.442 


100% 




0.6 


0.429 


97% 




1.5 


0.287 


65% 




2.3 


0.039 


9% 


4.95 





0.574 


100% 




0.6 


0.560 


98% 




1.5 


0.395 


69% 




2.3 


0.079 


14% 



TABLE II. T(Mf ax = 1.86GeV) with charged-lepton energy cuts. 



m b (GeV) Ef ut (GeV) T(M^ ax = 1.86 GeV)(|Kfe| 2 GeV) percentage 

4.8 0.442 88% 

0.6 0.432 86% 

1.5 0.312 62% 

2.3 0.054 11% 

4.65 0.359 82% 

0.6 0.351 79% 

1.5 0.254 58% 

2.3 0.039 9% 

4.95 0.538 94% 

0.6 0.526 92% 

1.5 0.381 67% 

2.3 0.079 14% 



TABLE III. T(Mf ax = 1.50GeV) with charged-lepton energy cuts. 

m ft (GeV) Ef ut (GeV) T(M^ ax = 1.50 GeV)(|Kb| 2 GeV) percentage 

4.8 0.348 69% 

0.6 0.342 67% 

1.5 0.261 52% 

4.65 0.263 60% 

0.6 0.259 59% 

1.5 0.200 45% 

4.95 0.465 81% 

0.6 0.457 80% 

1.5 0.343 60% 



TABLE IV. r(Mf ax ) with charged-lepton energy cuts (m b = 4.8 GeV and \i\ = 0.6 GeV 2 ). 



Mf ax (GeV) £f ut (GeV) T(M% ax )(\V ub \ 2 GeV) percentage 



no cut 0.519 100% 

0.6 0.505 97% 

1.5 0.351 68% 

1.86 0.460 89% 

0.6 0.450 87% 

1.5 0.329 63% 

1.5 0.391 75% 

0.6 0.384 74% 

1.5 0.290 56% 



TABLE V. r(M^ ax ) with charged-lepton energy cuts (m b = 4.8 GeV and a s = 0.35). 

M™ ax (GeV) ^ cut (GeV) r(M£ ax )(|K&| 2 GeV) percentage 

no cut 0.441 100% 

0.6 0.429 97% 

1.5 0.292 66% 

1.86 0.367 83% 

0.6 0.360 82% 

1.5 0.264 60% 

1.5 0.272 62% 

0.6 0.269 61% 

1.5 0.211 48% 



TABLE VI. r tota i with charged-lepton energy cuts (ACCMM distribution function of Eq. (21)). 



m 6 (GeV) Ef ut (GeV) r total (|Kb| 2 GeV) percentage 



4.8 0.503 100% 

0.6 0.489 97% 

1.5 0.336 67% 

2.3 0.055 11% 

4.65 0.443 100% 

0.6 0.429 97% 

1.5 0.287 65% 

2.3 0.040 9% 

4.95 0.573 100% 

0.6 0.558 97% 

1.5 0.394 69% 

2.3 0.079 14% 



TABLE VII. r(Mf ax = 1.50GeV) with charged-lepton energy cuts (ACCMM). 

m 6 (GeV) Ef ut (GeV) T{mT x = 1-50 GeV)(|Kfe| 2 GeV) percentage 

4.8 0.351 70% 

0.6 0.345 69% 

1.5 0.263 52% 

4.65 0.271 61% 

0.6 0.267 60% 

1.5 0.205 46% 

4.95 0.465 81% 

0.6 0.456 80% 

1.5 0.343 60% 



FIGURES 

FIG. 1. The 3-dimensional plot of the double differential decay rate dT / dEidMx with 
mj = 4.8 GeV. The lines in the bottom plane show the contour lines. The three straight lines in 
contour figure, which are the projections of three curved lines in the 3-dimensional figure, indicate the 
charged-lepton energy cuts, i.e. Ef ut = 0.6, 1.5,2.3 GeV, respectively. 

FIG. 2. Hadronic invariant mass distribution with m b = 4.8, 4.65, 4.95 GeV and Ef ut = 0, 0.6, 
1.5, 2.3 GeV. The four curves with the same line style are correspondsing to Ef ut = 0, 0.6, 1.5 and 
2.3 GeV, respectively, from top to bottom. 

FIG. 3. r(Mf ax ) with m b = 4.8, 4.65, 4.95 GeV and Ef ut = 0, 0.6, 1.5, 2.3 GeV, integrated up 
to Mjp". The four curves with the same line style are correspondsing to Ef nt = 0, 0.6, 1.5 and 2.3 
GeV, respectively, from top to bottom. 
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